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EFFtCT OF SEALS ON ROTOR SYSTEMS 


David P. Fleming 
NASA Lewis Research Center 
Cleveland, Ohio 44135 


Seals can exert large Torces on rotors. As an example, in turbopui^ 
ring seals film stiffness as high as 90 MN/m (500 000 Ib/in) have been 
calculated. This stiffness is comparable to the stiffness of rotor 
support bearings; thus seals can play an important part in supporting 
and stabilizing rotor systems. This paper reviews the work that has 
been done to dete»‘mine forces generated in ring seals. Working 
formulas are presented for seal stiffness and damping, and geometries 
to maximiz** stiffness are discussed. An example is described where a 
change in seal design stabil iz ed a previously unstable rotor . 


INTRODUCTION 


SYMBOLS 


Ring seals have received considerable study A 
over the past dozen years. The principal reason 
tor this attention is the Increasing importance a 

of seals in affecting the dynamic response of 
machinery. Large forces can be generated B 

between shafts and seals; these forces can be 
useful in supporting and stabilizing rotating 
machinery. Conversely, Improperly designed B 

seals can pronHjte instability. 


defined by Eq. {??) 
defined by Eq. (?3) 
seal damping coefficient 

BC 

dimensionless damping, ->i 


Figure 1 lllustreates a ring seal. This 
seal has the appearance of a journal bearing, so B 
it Is not surprising that bearing- like torces 
can be generated. Seals typically have larger 
clearances than bearings; thus hydrodynamical ly .. 

generated forces would be lower. (Stiffness, 
which equals tone/ displacement, changes as the 
Inverse cube of clearance for laminar flow. , 

Stiffness is frequently used rather than force ^ 

to characterize a seal or bearing.) However, 
the pressure drop across a seal yields a 
hydrostatic force which varies directly with the 
pressure drop. This force can be quite large in 
high-pressure machliiery; for example, the ^ 

intc’'.tage seals of the Space Shuttle fuel pump, 
with a pressure drop of 140 bars (i’UOO psi) has F 

a calculated stittness of iS MN/m (,’.'0 000 
lb/ in.) [1]. G 


dimensionless damping, -k 

L D 

radial clearance of concentric seal 


dimensionless damping, B 



seal diameter 
defined by Eq. (21 ) 
seal radial force 
seal inert la coet t iclent 


The purpose of this paper Is to review the 
developim'iit of force determination In ring 
seats, and to present the latest, most accurate 
methods of calculating seal torces. 
Additionally, seal geon«’tries to maximize 
generated torces will be discussed. An example 
of seal effect on rotordynamic behavior will 
also be presented. 


h local seal clearance 


k seal stiffness 


K dimensionless stiffness, — qj 

'’o 






L seal length 

M leakage flow 

p pressure 

q Rw/U 

R sci’> radius 

Re Reynolds number, 2CpU/v 

Re^ "sonic" Reynolds number, 2C 

u local fluid velocity 

U u for centered seal 

V swirl velocity 

H swirl ratio V/Ru 

x,y transverse coordinates 

z axial coordinate 

a seal taper half-angle 

Y specific heat ratio 

c eccentricity ratio e/c 

n total entran:e loss factor 1 + t 
X friction factor 

M fluid viscosity 

i entrance loss factor 

p fluid density 

o dimensionless seal length x!./C 

u) angular velocity 

Subscripts: 

0 upstream stagnation region 

1 seal entrance 

2 seal exit 

ANALYTICAL PRINCIPLES 

In the earliest analyses numerous 
simplifying assumptions were made. Later work 
removed some of these assumptions. In order to 
aid in understanding how a high-pressure seal 
operates, it is useful to go through the most 
elementary analysis. For this analysis the list 
of assumptions is, of course, the longest and 
reads 


1. The sealed fluid Is incompressible. 

2. Effects of rotation are neglected. 

3. Fluid flow Is one-dimensional In the 
axial direction, that Is, any circumferential 
flow Is neglected. 

4. The eccentricity is small compared with 
the seal concentric clearance. 

5. The fluid friction factor Is constant 
throughout the seal. 

6. All flow is steady state. 

7. Entrance losses can be accounted for by 
a suitable choice of loss coefficient. 

8. There is no pressure recovery at the 
seal exit. 

The analysis begins with the continuity and 
momentum equations. These are given in [2] and 


for invariance with time are 

h ^ = 0 (1) 

The boundary conditions are 

P = Pq - npu^/2 at z = 0 (3) 

p = 0 at z = L (4) 


where n is the total entrance loss factor. It 
is che sum of the pressure drop, pu 2/2, that 
occurs due to Bernoulli's principle, and the 
pressure drop due to the developing velocity 
profile in the fluid. The value of n is 1.65 
for laminar flow [3] and drops to near 1 for 
highly turbulent flow [4] 

The solution to Eq. (1) is 

u = constant (5) 

and to Eq. (2) 



Applying the boundary conditions, 

Po " r) 

P •= Pq(L - (8) 

The total restoring force on the seal can be 
found by integrating the restoring component of 
pressure over the seal area 
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rl 

F - -R / I p cos « d 2 (to (9) 

*'0 •'0 

The local clearance h Is given by 

h - C(1 ♦ c cos e) (10) 

Invoking the small eccentricity assumption 
(c « 1), and substituting Eqs. 

(8) and (lO) In (9), 


pressure drop and again the pressure profile 
changes little with eccentricity. Thus the 
greatest stiffness occurs when the pressure drop 
due to entrance effects Is comparable to that 
due to Mall friction. This Is similar to the 
situation In externally pressurized bearings, 
Mhere for maximum stiffness the restrictor 
pressure drop equals the film pressure drop. To 
carry the analogy further, one could say that 
the seal entrance Is analogous to the restrictor 
In a pressurized bearing. 


F - w*Rp^L^nC*/(nC ♦ 2xL)^ 


( 11 ) 


or, defining a dimensionless stiffness 






( 12 ) 


Mhere 


o • xL/C 


(13) 


A SHORT HISTORY OF ANALYTICAL DEVELOPMENT 

The foundations of seal force analysis Mere 
laid by Henry Black and coMorkers [5,6,7], In 
keeping with the assumptions listed above. 
Black's analyses are one-dimensional. 
Essentially, Black looks at flow In a duct Mith 
a wall which moves transversely; thus the flow 
Is unsteady. The result; consist of 
coefficients for stiffness, damping, and 
Inertia, such that the force on the seal Is 
given by 

F - kx + Bi + Gx (14) 


Having derived this simple expression for seal 
stiffness, we can use It and the seal model to 
gain some understanding of how the sea’ develops 
its stiffness. Figure 2 shows the fluid 
pressure In the concentric seal as a function of 
axial distance. After an Initial abrupt drop 
(due to the Bernoulli effect and entrance loss, 
Eq. (3), the pressure drops linearly to zero at 
the exit. The clearance and fluid velocity are 
uniform around the circumference. 

Now consider what happens when the seal Is 
moved to the side. Fig. 3. The clearance varies 
around the circumference according to Eq. (10). 

On the high-clearance side, the fluid can move 
faster than before, and, conversely, on the 
low-clearance side If has to go slower. The 
entrance pressure drop depends on the fluid 
velocity, so there's a greater Initial drop on 
the hlqh-c learance side. On both sides, the 
pressure must drop to zero at the seal exit. 

Thus the pressure profiles look like those 
depicted In Fig. 3. The pressure Is general 1» 
higher on the low-c learance side, and this 
results in a force that tries to center tne seal. 

Figure 4 Is a plot of the dimensionless 
seal stiffness, k, versus the parameter o 
which may be thought of as a dimensionless seal 
length. Ihe entrance loss factor was chosen as 
ii - 1.1. wtii(h IS a value representat Ive of 
turbulent flow. For low values of a (i.e., 
short seals) stiffness rises with seal length. 

A peak IS reached, and then stiffness drops as 
length Increases. This phenonwion may be 
explained as follows: for very short seals, the 

entrance effect is the domtnant pressure drop 
mechanism, that is, the seal behaves somewhat as 
an orifice. With wall friction playing such a 
small part , the pressure profile does not change 
much with seal eccentricity (see Fig. 5) and 
thus the stiffness is low. Conversely, for long 
seals wall friction causes the bulk of tfie 


In [6], shaft rotation is accounted for by 
assuming that the entire flow field rotates at 
one-half the shaft speed. The fluid equations 
are solved in a coordinate system which also 
rotates at half shaft speed (with fluid being 
assumed to flow axially in the rotating 
system). A coordinate transformation then 
yields an expression for forces on the shaft In 
the X and v directions 



B 

-Gw 



K, B, and G are the same as those for the 
nonrotafing case. Black did, however, modify 
the friction factor to account for tne 
additional fluid velocity due to rotation. For 
turbulent flow 



This rn>"m was suggested by Yamada [8]. It may 
be noted chat the nx'an flow velocity U Is not 
known ajgri^orl; It must be found through 
IteratTve scTTut ions of Eqs. (7) and (lb). In 
[7], the c ircumterent ial variation of trictlon 
factor Is accounted for; this results in a 
substantial increase in ttie direct stiftness 
coefficient kxx 'd- b) but relatively small 
changes in the damping and inertia coet t ic ients , 
compared with the constant trictlon factor 
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solutior.. As Mg. b shows, the rotational 
velocity affects the results; stiffness rises as 
rotational speed increases and, as shown in [7], 
damping decreases. 

The works cited up to now have all used a 
small eccentricity analysis and presented 
results for small perturbations about a centered 
seal position. Allaire, Lee, and Gunter [2] 
extended the analysis to obtain results for a 
perturbation about a finite eccentricity 
position. They employed numer'cal integration 
to obtain stiffness, damping, and total load 
capacity for various seal eccentricities. In 
another departure from earlier analyses, the 
authors of [2] neglected the time variation of 
fluid velocity in the governing equations. This 
change simplifies the solution procedure but 
inevitably introduces some inaccuracy in the 
calculated damping coefficient (the stiffness is 
unaffected). In [Ij this inaccuracy wa-" ’ .id 
to be no more than 16 percent, with th« 
simplified solution overestimating the oamping. 

The results of the finite eccentricity 
calculations showed that, in general, seal 
stiffness and damping are not much affected by 
eccentricity (Figs. 1 and 8). The exception is 
the damping in the direction of displacement, 
where, as would be expected, the oamping 
coefficient increases considerahlv as the seal 
journal approaches the side wall. Seal rotation 
was not considered in [?J. 

Effec t of inlet swirl. - As noted above, 
Black Fb] acc'o’uVi t'ed' T or sea 1 rotation by 
assuming that the fluid circumferential velocity 
is one-halt of the shaft speed throughout the 
seal. Actually, the fluid velocity only 
approaches this limit asyn^tot ical 1y as it 
proceeds ti'-ouoS the seal. In applications 
such as centritugal pum;' interstage seals, the 
fluid may enter the seal with little or no 
rotational velocity. Furthenm're, it the seal 
length is short, thiid rotational speed may 
differ appreciably from the halt-shaft-speed 
asymptotic value even at the seal exit. Thus 
seal dynamic properties are likely to vary with 
the amount ot inlet swirl present. The effect 
of inlet Swirl has been investigated in two 
papers, one by Black, Allaire, and Barrett [9] 
and one by Childs (10]. The analysis in [10] is 
the more comprehensive, as it solves the 
short-bearing turbulent lubrication equations ot 
Hirs [11,12]. Since several ot the assumptions 
listed in the Analysis section need no longer be 
made, it is woFlTiwIilTe to list those assumptions 
still required by [10]. 

1. The sealed fluid is inconniressible. 

Pressure induced circumferential flow is 
neglected (the "short bearing" app'^ox in.af ion 

[i-i]). 

3. The eccentricity is small cimipared with 
the seal concentric clearance. 


4. Entrance losses can be accounted for by 
a loss coefficient. 

5. There Is no pressure recovery at the 
seal exit. 

£h11ds results [10] for direct and cross-coupled 
stiffness and damping are given for turbulent 



Also, q is defined as the ratio of journal 
rotational speed to mean axial velocity. 


The initial swirl ratio Wj is the ratio Ot 
circumferential fluid velocity at the seal 
entrance to the journal spred. 
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The swirl ratio approaches 1/2 as the fluid 
proceeds through the seal; Ui-1/2 implies the 
half-journal speed assumption of Black and 
Jenssen [6]. 

Whit. e expressions above appear 
comp Heater they can be programmed on a pocket 
calcula;’'- ' digital computer in a 
straiohtti.rward manner. Figure 9 shows the 
results of the calculations. The seal used for 
these calculations has a length to diameter 
ratio of 1/2; the circumferential journal speed 
is 1/2 that of the mean axial fluid velocity. 

As expected, prnswirl has a stronger effect 
for small values of a, the dimensionless seal 
length. The larger the initial value of 
circumferential fluid velocity, the smaller is 
the direct stiffness and the larger the 
cross-coupled stiffness and damping, while 
direct damping is completely unaffected. 
Cross-coupled stiffness is usually viewed as 
promoting instability. Small values of preswirl 
reduce l^xy and thus would appear to be 
stabilizing. The authors of [9] point out that 
negative preswirl - the fluid revolving opposite 
the direction of shaft motion - can result in 
negative Kxy and thus counteract 
whirl-producing forces elsewhere on the shaft. 

State of the art . - This term is used with 
some misgivings, as ring seals continue to 
received intensive study and new publications 
appear frequently. 

Equations (17) to (20) presented above 
appear to be the most comprei'',°nsive available to 
calculate seal stiffness and a.mping, requiring 
the fewest simplifying assumptions and 
presetting results in explicit form. Results 
are imited ti small displacements from the 
centered posiMon, but, as shown in [2], seal 
stiffness and damping are largely insensitive to 
seal eccentricity. The source of Eqs. (17) to 
(20), [10], also presents expressions for 
inertia coefficients, but these are usually of 
minor importance. 

Recently two authors have investigated 
seals of finite length (thereby eliminating the 
short seal assumption). In [K] the turbulent 
Reynolds equation was solved numerically for 
finite eccentricity ratios, using the finite 
element method. Childs, in [15], starts with 
Hirs' bulk flow equations [11,12], transforms 
them to orditiary differential equations by small 
eccentricity assumptions, and integrates 
numerically. Both references report that the 
short seal solution overestimates the seal 
coefficients with the error increasing for 
longer (L/P) Seals. 

In this connection, Black and Jenssen [b] 
P'"oposed correction factors for finite length 
seals based on an approximate analysis. The 
results of [14] and [15] indicate that these 
factors result in overcorrft ion . A nwre 
accurate formula would seem tc. he 



EXPERIMENTAL CORRELATION 

Experimental data, particularly dynamic 
data, are very sparse. Black and Jenssen show 
results of experiments in [6] and [7] for axial 
Reynolds numbers up to 20 000 and ro' jtional 
Reynolds numbers (based on half journal speed) 
up to 14 000. Both unidirectional pnd rotating 
load (due to unbalance) tests were run. In both 
papers, agreement between analysis and 
experiment is shown as close, although the 
analysis of [7] employed the local friction 
factor effect while that of [6] did not. Childs 
and Dressman [16] have presented preliminary 
results which show that measured direct 
stiffness is generally a good deal higher than 
predicted by [10] (to 90 percent) while direct 
damping is much lower than predicted. 

OPTIMUM GEOMETRIES 

Thus far it has been tacitly assumed that 
the seal clearance is constant in the axial 
direction; in other words the seal has a 
straight bore. From physical reasoning it is 
easy to see that the direct stiffness would be 
higher if the clearance at the fluid exit were 
less than at the entrance. In [17] Fleming 
calculated stiffness for seals with stepped and 
tapered bores (Fig. 10) and determined optimum 
geometries to maximize either the stiffness K 
on the ratio of stiffness to leakage K/M. This 
was followed with calculations for dancing in 
tapered bore seals [1]. Figure 11 shows that 
substantial increases in stiffness are possible 
with an optimum taper configuration. The 
stiffness increase is obtained with only a 
moderate leakage penalty. The inlet-to-out let 
clearance ratio for the maximum K/M seal is 
nearly constant over the entire range of 
dimensionless seal lengths (o) shown in Fig. 11 
having a mean value of 1.9. Direct damping is 
lower in tapered seals as shown in Fig. 12. 

This is principally because of tne tapered seals 
larger average clearance (tapered and straight 
bore seals are compared using the minimum 
clearance in the seal). More recently Childs 
extended his solution of Hirs' lubrication 
equations to tapered bore seals [18]. His 
results agree qualitatively with those of [1] 
and [17], Childs considers variable inlet swirl 
and rotational effects within the constraint of 
the short bear ing appi'Ox iiiiat uin; thus greater 
accuracy would be expected. 

SEALS FOR COMPRESSIBLE FLUIDS 

The flow of a compressible fluid (most 
commonly a gas) differs from that of a-, 
incompressible fluid in that change:, in pressure 
are accompanied by changes in pressure and 
temperature. For flow through a seal - usually 
approximated as adiabatic - as the p- essure 
drops along .'he flow path the dens'ty and 
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tMperature «Uo drop, and the velocity 
Increases. Except for the special case of a 
convergtn^diverglng channel, the maximum fluid 
speed 1$ limited by the speed of sound. Thus 
for some ratio of upstream to downstream 
pressure the flow will become choked and hence 
unaffected by further decreases In downstream 
pressure. 

Fluid compressibility has a large effect on 
the forces generated In seals. Stiffness and 
damping coefficients are calculated for both 
straight bcre and tapered-bore seals In [ 4 ] and 
[19]. All of the assumptions listed In the 
analysis section apply, except, of course, that 
compressibility Is allowed; the fluid Is assumed 
to be a perfect gas. Figure 13 shows stiffness 
of straight seals and optimum tapered seals for 
choked flow. The abscissa Is the seal 
clearance-to-length ratio C 2 /L times a 
Reynolds number Re» formed using the speed 
of sound for conditions upstream of the seal. 

RCq - 2C^^oPo/uo 

The curves on the left are for laminar flow. 

The nearly horizontal curves on the right are 
for turbulent flow; there Is a separate curve 
for each value of C 2 /L. The left end of the 
curves for turbulent flow corresponds to a 
Reynolds number In the seal passage of 3000; 
this Is generally considered the lowest value 
for which one can be assured of turbulent flow. 

A Reynolds number of 2300 1s usually taken as 
the upper limit for laminar flow. Points where 
Re - 2300 are shown for various C 2 /L values 
on the laminar flow curves. 

The surprising feature of the results is 
the negative stiffness predicted for low values 
of C 2 /L. This is analogous to the lockup 
phenomenon observed in pressurized 
gas-lubricated bearings. Negative stiffness 
could cause rapid failure of floating ring seals 
anJ would rarely be beneficial even for rigidly 
mounted seals. Tapered seal stiffness is always 
positive and generally much higher than straight 
seal stiffness. The optimum clearance ratio is 
near 1.9 for most conditions, as it is for 
incompressible fluids. 

Damping is shown in Fig. 14. As for 
stiffness, laminar flow curves are on the left 
and turbulent-flow curves (a separate one for 
each C 2 /L) are on the right. Again, as for 
incompressible fluids, damping is lower for 
tapered seals. Damping is minimized for 
c Isirancc-tc- length ratios of 0. 002-0, OOS and is 
higher for ratios outside this range. 

Expermental data are again sparse. 

Hendricks [20] obtained pressure measurements 
along the length of straight and tapered bore 
seals in the concentric and fully eccentric 
positions. Because of the small number of 
pressure taps he did not attempt to determine a 
net forte generated by the seal; however, he 
estimated that the overall stiffness of the 


tapered seal was no higher than that of the 
straight seal. Physical reasoning leads one to 
question this conclusion. 

In work by Burcham and Diamond [21] the 
tapered bore seal was far superior to an 
alternative seal, a Rayleigh lift pad design. 
Both were floating ring seals, and the objective 
was to determine the endurance of the seals 
while maintaining satisfactory leakage control 
(I.e., minimizing wear). No Internal pressure 
measurements or force measurements were made. 

In this type of test, the seal having a higher 
stiffness (higher centering force) can better 
minimize the rubs which cause wear. The 
Investigators noted that there was less seal 
wear when sealed pressure was higher. This is 
In agreement with the analysis, which shows that 
a higher sealed pressure produces higher sea' 
stiffness. 

SEAL EFFECT ON ROTOR DYNAMICS 

When significant forces (or equivalently, 
stiffnesses) are generated by seals, the effect 
on the shaft-rotor system is the same as if a 
bearing of that stiffness were placed on the 
shaft. Consequently, changes in dynamic 
response can be expected as a function of seal 
stiffness. 

An example of seal forces being exploited 
to favorably influence rotor dynamic response is 
the hign pressure fuel pump for the space 
shuttle main engine. The rotor of this pump is 
shown schematically in Fig. 15. Three 
centrifugal pump stages are driven by a 
two-stage axial turbine. The rotor is supported 
on two duplex ball bearings. Ring seals are 
used between pump stages. Originally, these 
seals were of serrated design (Fig. 16) to 
deliberately reduce the forces produced, as it 
was believed the cross-coupled stiffness would 
promote instability. In fact, other exciting 
foices produced a subsynch-onous instability 
which tegan as the pump speed passed twice the 
first critical speed. Various modifications 
were tried to eliminate the instability. There 
was little success until the serrated seals were 
replaced with snwoth-bore seals [22]. With this 
change the pump could be run to full speed 
without exhibiting signs of instability. The 
conclusion drawn was that the higher direct 
stiffness of the smooth seals was largely 
responsible for the stabi 1 izatioti and that the 
benefit of higher direct stiffness overshadowed 
the deleterous effect of the higher 
cross-coLinled stiffness. The effect of 
higher-stiffness seals can be quantified in one 
respect by determining the effect on 
resonances. This is shown in Fig. 17 where 
resonant frequency is plotted as a function of 
interstage seal stiffness tor a rotative speed 
of 37 000 rpm, the normal operating speed. The 
stiffness of the original serrated seals is not 
known precisely; however seal stiffness has 
little effect on resonant frequency for low 
values of stiffness. The first and second 
resonances are 9500 and 19 8tX) cpm. 
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respectively. When the serrated seals are 
replaced by smooth-bore seals, the resonances 
have risen to 13 $00 and 23 000 cpm. Further 
Increases occur If optimum tapered-bore seals, 
described earlier, are used. The two resonances 
then are 16 000 and 28 000 cpm. These 
frequencies are 70 and 40 percent higher, 
respectively, than the resonant frequencies with 
the original serrated seals. 

CONCLUDING REMARKS 

Significant forces can be generated by ring 
seals; film stiffnesses as high as 90 MN/m 
(500000 1b/ In.) have been calculated for 
turbopump ring seals. These forces have been 
effectively utilized to support and stabilize 
shaft-rotor systems. This paper has presented a 
summary of work that has been done to determine 
forces generated In ring seals. Both 
compressible and Incompressible fluids were 
considered; working formulas were presented to 
calculate seal stiffness and damping. Seal 
geometries to maximize stiffness were 
described. An example showed the effect seals 
can have on the dynamic behavior of rotor 
systems; use of high-stiffness seals raised the 
rotor critical speeds by nearly 50 percent. 
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Fiqur* 3. - Velocity and preiiure profile in eccentric 
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Figure 9. - Continued. 
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Figure 9. - Concluded. 











